NASA TN D-3 


loan COPY: RFTUR,E = 
AFWL fWLlL-2] 

K,Ki land afb, n □= 


A THEORY FOR INFLATED 
THIN- WALL CYLINDRICAL BEAMS 


by W. B. Fichter 

Langley Research Center 
Langley Station, Hampton, V a. 



tK A S '°V 

// V 
'a/.;/;... ' 

^ va 


V ^.*.7 


s*sZ 

WASHINGTON, 0. C. . JUNE 1966 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


IRARY KAFB, NM 




TECH LIBRARY KAFB, NM 



□13DM1D 

NASA TN D-3466 


A THEORY FOR INFLATED THIN-WALL CYLINDRICAL BEAMS 

By W. B. Fichte r 

Langley Research Center 
Langley Station, Hampton, Va. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 — Price $1.00 


A THEORY FOR INFLATED THIN -WALL CYLINDRICAL BEAMS 


By W. B. Fichter 
Langley Research Center 

SUMMARY 

Nonlinear equilibrium equations are derived for the twisting, bending, and 
stretching of pressurized thin-wall cylindrical beams. Thin-wall beams, including those 
which depend on internal pressure for their load-carrying ability, appear to have appli- 
cation where lightweight structural members are required. Example problems involving 
column buckling and the bending of a beam column are solved with a set of linearized 
equations to illustrate the theory. The linearized equations reduce to the Timoshenko 
beam equations when internal pressure and axial force are set equal to zero. 

INTRODUCTION 

Lightweight inflatable plates and cylinders which depend on internal pressure for 
much of their load-carrying ability are receiving considerable attention in aerospace 
research. Packageability and automatic erectability make inflatable construction espe- 
cially attractive for satellites and space stations. Inflatable construction may also pro- 
vide a suitable means of building reentry vehicles with wing loadings low enough to 
alleviate significantly the problem of reentry heating. In addition, inflatable membrane 
cylinders may be quite useful as booms in a paraglider or similar vehicle. 

For all these applications, information is needed on the load-carrying capabilities 
of inflatable structures. Simple methods of analysis have been used successfully to 
determine collapse loads for certain types of inflated structures (refs. 1 and 2). In addi- 
tion, a linear theory for the structural analysis of inflatable plates has been presented in 
reference 3 and has been compared favorably with experiments in reference 4. Refer- 
ence 5 contains a combined aerodynamic and structural analysis of a paraglider wing 
consisting of a flexible triangular sail between three equally spaced booms which are 
joined at the nose. In the analysis, stresses and deflections of the sail were calculated 
under the assumption that the booms were rigid. The results of the present paper may 
be of use in extending the analysis of the paraglider wing with rigid booms to the case of 
a wing with flexible booms. 

The purpose of the present investigation is to develop for the long, inflated, 
circular -cylindrical beam a structural theory which allows for moderately large 



displacements, and which can be reduced in a consistent manner to a simple system of 
linear equations. The theory is applicable to pressurized membrane cylinders (that is, 
cylinders having negligible wall bending stiffness) when the internal pressure and applied 
forces are such that the cylinder wall remains in tension. The application, however, is 
not limited to membrane cylinders; the theory can be applied to pressurized or unpres- 
surized thin -wall cylindrical beams in which wall bending stiffness is sufficient to pre- 
vent local buckling of the wall. 


SYMBOLS 

A cross-sectional area of cylindrical shell, 27rrt 

BpBg arbitrary constants 

C constant of integration 

E Young's modulus 

F compressive axial force 

G shear modulus 

I area moment of inertia of cylindrical shell, irr^t 

2 

Kxa Euler buckling load, — 77 

a 9 

Kg shear stiffness parameter, G 7 rrt 

L length of cylindrical beam 

MpMgjMg moments about the X-, Y-, and Z-axes, respectively (see fig. 1) 
N total axial force on a beam cross section (see fig. 1) 

N x ,N xe axial- and shear-stress resultants in shell 
P axial force due to internal pressure, pwr^ 
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p 


internal pressure 


Q payload 

applied distributed loads in the X-, Y-, and Z-directions, respectively 
r radius of circular cylindrical shell 

T^Tg,!^ applied distributed moments about X-, Y-, and Z-axes, respectively 

t cylinder wall thickness 

translations of abeam cross section in X-, Y-, and Z-directions, 
respectively 

u,v,w axial, tangential, and radial displacements of a point in surface of shell 

(see fig. 1) 

VgjVg shear forces in Y- and Z-directions, respectively (see fig. 1) 

AV change in enclosed volume of cylindrical shell due to deformation 

W work done by externally applied loads during deformation 

X,Y,Z rectangular Cartesian coordinates (see fig. 1) 

x,y,z distances along X-, Y-, and Z-axes, respectively 

0 parameter defined in problem 1, /3 2 = — — — + ( ^ 7rrt — 

ETrr^tVP + G7rrt - F 

y x9 local shear strain in membrane 

7 1 angle of twist of beam cross section per unit length 

y 2>73 shear strain of abeam cross section in the Y- and Z-directions, respectively 

e x local extensional strain in shell in X-direction 
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e ^ extensional strain of cylindrical beam in X-direction 

£, 77 , £ coordinates of a point in deformed cylindrical shell, referred to the X-, Y-, 

and Z-axes, respectively 

6 circumferential angular coordinate 

k.Q ', Kg curvature of cross-section center line in XZ- and XY -planes, respectively 

„ N(P + K s ) 

X parameter defined in problem 2, X = 

EI(P + K s - N) 

IIj strain energy of deformed cylindrical shell 

Ilg change in potential energy of pressurizing gas due to deformation 

Pq,Pl variables of integration (see eq. (7)) 

o>i,o> 2 ,u >2 rotation of a beam cross section in the YZ-, XZ-, XY-planes, respectively 
Subscripts: 
cr critical 

i index 

Barred quantities represent the applied forces or moments at the ends of the beam. 

ANALYSIS 

A sketch of the inflatable circular -cylindrical beam under investigation, showing the 
forces, moments, displacements, and coordinate systems, is presented in figure 1. The 
ends of the beam are considered to serve no function other than to contain the internal 
pressure. 

Two important simplifying assumptions are made in the analysis: first, any cross 
section of the pressurized cylindrical beam remains undeformed under a system of loads, 
so that deformations of the beam are restricted to translations and rotations of a beam 
cross section; and second, the translations and rotations of a cross section are small 
enough that the displacements of a point in the surface of the beam can be represented 
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with little error by the vector sums of the displacements due to translations and rota- 
tions of a cross section. A further simplification attendant to these assumptions is that 
the circumferential strain is negligible and may be disregarded in the analysis. Under 
the present assumptions, the theory holds for long cylindrical beams having walls with 
either significant or negligible bending stiffness, provided that local buckling or wrinkling 
does not occur. 


Derivation of Governing Nonlinear Equations 
In order to derive a system of equations, use is made of the variational equation 

6(n 1 + n 2 - W) = 0 (1) 

where II ^ is the strain energy of the cylinder, Ilg is the change in the potential energy 
of the pressurizing gas due to deformation, and W is the work done by externally 
applied loads during deformation. If the circumferential strain is neglected, then 


•i pL /*» 277 

n i“jj 0 Jo ( n * 6 x + N xeW) r de * 


( 2 ) 


Also, 


n 2 = - P AV 


(3) 


where p is the internal pressure (assumed constant during deformation) and AV is 
the change in enclosed volume due to deformation. 

From reference 6, the nonlinear strain-displacement relations which include all 
second-degree terms are found, after some manipulation, to be 


e = — + - 
x 9x 2 


8V \ 2 + /0w\ 
9X ) \9X j 


(4a) 


and 


_ 9v_ 9u\ 1_ 9U/. _ 1 9v w\ 1 9w/9w \ 

7 x0 9x \ ~ 9x / r 9 0 \ r 90 " r J r 9x \90 ' / 


However, with the assumption that 
shear strain reduces to 



are small compared with 1, the 
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r -9V + I9H+19W3W x 

7 x0 9 X r 90 r 9x \80 ) 

In terms of the translations and small rotations of a cross section, the local displace- 
ments are given by: 

u = Uj + r sin o>g cos 0 - r sin o>g sin 0 

v = Ug cos 0 - Ug sin 0 - r sin + |-(cos Wg - cos cog^sin 20 
w = U 2 sin 0 + Ug cos 0 - r jl - cos + (1 - cos Wg^cos^ + (1 - cos Wg^sin^Jj 
For small rotations 


sin Wj ~ 


cos co. 25 1 — 


and the local displacements reduce to 


u = U x + r u >2 cos 0 - rcog sin 0 


v = U 2 cos 0 - Ug sin 0 - rco^ + --^ajg^ “ co 3 ^) s ^ n 2< 


w = Ug sin 0 + Ug cos 0 - + cog^cos^ + Wg^sin^^ j 

The change in enclosed volume due to these displacements may now be computed by 
expressing the coordinates of a point in the deformed cylindrical beam in terms of the 
local displacements and the coordinates of the undeformed beam. Let 


| = u + x 


77 = (r + w)sin 0 + v cos 


£ = (r + w)cos 0 - v sin 0 


Then by the use of Gauss' theorem (ref. 7) the change in volume is 


S L «“> 


3(x,0) 


dx d0 
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where the integration is over the total deformed surface S. By a transformation of 
coordinates, 


AV = -tit h + 


f 2 * f L 77 ) 4^4 dd + f 27F r dp L d0 

J 0 J 0 9(x,0) J 0 ^0 ^ L L ^(P l , 0 l ) L 


p27r p r 

S 0 (h%) - \ dp o de Q 

J 0 Jo ^ u ° ; 9(P O ,0 O 0 ° 


( 7 ) 


where the subscripts 0 and L refer to the ends x = 0 and x = L, respectively. 

Inasmuch as the cross-section displacements and rotations are functions only 
of x, integration with respect to 9 and substitution from equations (5) and (6) yield 

.L 


AV = ”' 2 I 0 [< + U2 '“3 ' U 3 “2 ' 2 (" 2 2 + “ S 2 )]* 5 


(8) 


where primes denote total differentiation with respect to x, and where terms of the third 
or higher degree in the displacements and rotations and their derivatives are neglected. 

To the same degree of approximation, the local strains expressed in terms of 
cross-section rotations and translations are obtained by substituting equations (5) into 
equations (4): 


e x = Uj' + rw 2 'cos 9 - raisin 9 + 


- 2r 


fEv 

'cos 9 - sin 




2 *2 o * 2 

+ Ug + r^a^ 


/ ( 9 ) 


?x9 


= -ro> 


1* + ( U 2* - w 3 + W 1 U 3’) 


OJg + )COS 6 

co 0 - o) 0 f (U 0 ^sin 26 


( u 3 ’ + w 2 - w^’jsin 9 


K-2-2 - w 3' w 3) si 


Substitution of equations (9) into equation (2) gives 
>2 tt 


n i = I I 0 ^ ( N x[ U l’ +r K - w l'u 2 ’)cos 9 - r(a) 3 ’ - c^'Ug’) 


sin 9 


IL2.. .2 + y^'2 + y^'2 


If 2 
2 V W 1 


)J + N x0 [- rco l' + ( U 2* - w 3 + w iU 3 ')cos 0 


- 3 ' + cl>2 - '^sin 9 + ^(^2 ^2 " a, 3' w 3] sin 26>|Jdx d0 (10) 
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In terms of the forces and moments on a beam cross section (see fig. 1), the 
stress resultants are defined by 


N ^2 Mq 

N„ = -1— + — =- cos 9 ^ sin 9 

irr 2 


x 2irr 7n*2 


M, V 9 vo 

N fl = rr + —cos 9 sin 6 

xB 2nr 2 


V 

-cos 9 

irr irr 

Substitution of equations (11) into equation (10) and integration over 6 yield 

n l = 2 I 0 ^ U l’ + + U 2 2 + U3 ’ 2 ) + + M 3 {^3 - "l’ 1 

+ MjWj’ + V 2 (U 2 ’ - + 0)^3*) + V 3 (U 3 ' + co 2 - WjUg'^dx 

If now direct, bending, and shearing strains for the beam are defined by 

tt ' 1 / 2 *2 t2 1 2 \ 

e l = U l + 2 j r “l + U 2 +U 3 j 
k 2 = w 2 ’ - 0)^2* 

“3 = u 3 ~ "l’^V 


then 


y i = w i 


^2 = U 2* - w 3 + "lV 


y 3 


= U, + co 0 - w.U, 
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n i = \ I 0 ( Ne l + M 2 x 2 + M 3 k 3 + M l*l + ¥2 + Vs)* 
In terms of the beam strains, the forces and moments may be defined by 


local 

(ID 

V) 

( 12 ) 

(13a) 
(13b) 
(13c) 
(13d) 
(13e) 
(13f ) 

(14) 
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N = 2E7rrte^ 

(15a) 

V 2 = G7rrty 2 

(15b) 

V 3 = Gffrtyg 

(15c) 

Mj = 2Girr\y^ 

(15d) 

3 

M 2 = Eirr t/c 2 

(15e) 

3 

Mg = Eirr t Kg 

(15f) 


Because of the assumption that the circumferential strain is zero, the material constants 
E and G may be defined independently and, hence, the definitions of these constants 
hold for orthotropic materials. Therefore, the variation of the strain energy may be 
written as 


OT 1 = In ( N6£ l + “2 6*2 + MgS/Cg + Mj6y^ + V 2 6y 2 + ^3^3 j dx 


(16) 


The variation of work due to externally applied loads may be defined by 
n L 

5W = j + + q 3 5U 3 + T i 6u> i + T 2 5ct> 2 + T 3 6 ct) 3] < ^ x 


( j 


+ ( N6Uj + + Vg6Ug + + M 2 5u> 2 + Mg6c0g^ 


(17) 


Substitution of equations (3), ( 8 ), (16), and (17) into equation (1) yields the variational 
equation 


>( n 1 + n 2 -w)=n(N & 1 + v 2 6 y 2 + Vg 6 y 3 + Mj 6 yj + M 2 6 k 2 


+ “s 6 ^) - P’ r26 [ U l' + U 2 '"3 - U 3“2 - j(“2 2 + “3 2 ) 

- + q 2 5U 2 + q 3 6U 3 + T^Sw^ + T 2 6u> 2 + Tg6u>gjjdx 

- (N6U 1 + V 2 5U 2 + V 3 6U 3 + MjScOj + M 2 6u> 2 + Mg6a>g\ L = 

f r\ 


( 18 ) 
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Use of the fundamental lemma of the calculus of variations, and integration by parts 
where necessary, yield the equilibrium equations and boundary conditions: 


N’ + q : = 0 


( 1TO 2')'-( M 2" l ')' + V 2'-( V 3“l') 

\ - Pa V + ^2 = 0 

(™3’)' -( M 3“l')' + V 3 , + ( V 2“l') 

’ + + ^3 = 0 

r^N^y + M,' -(M 2 U 2 ' + M 3 U 3 ') 

1 + V 3 U 2 " V 2 U 3 + T 1 = 0 

M 2'- V 3- P (“2 + U 3') + T 2 = 0 


M 3 , + V 2 + P ( U 2' -"s) + T 3 =0 

J 


where P = p7rr2. The boundary conditions at the ends of the beam are either 

0 or N - P - N = 0 

0 or NU 2 ’ - M 2 w 1 ’ + V 2 - V 3 w 1 - Pw 3 - V 2 = 0 
0 Or NUg’ - MgWj' + Vg + V 2 W 1 + PU> 2 - Vg = 0 



w 1 = 0 or r 2 Nw 1 ’ + M 1 - MgUg’ - MgUg' - Mj = 0 


(20) 


co 2 - 0 or M 2 - M 2 = 0 


ci >2 = 0 or Mg - Mg = 0 


Equations (19), along with equations (13) substituted into equations (15), constitute a 
system of 12 equations in the 12 unknown displacements, rotations, forces, and moments. 
As long as the nonlinear terms are retained, solutions usually will be obtained only by 
numerical integration or by one of several approximate -solution methods. 

Linearized Equations 

For many practical applications involving pressurized cylindrical beams, twist 
about the longitudinal axis is negligibly small. With the assumption of zero twist about 
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the X-axis, equations (19) and (20) simplify to 


N' + q x = 0 



(21a) 

( NU 2' j + V 2 ’ 

- Pcu 3 ’ + 

O 

li 

CM 

& 

(21b) 

( NU 3 ')’ + V 3' 

+ Pa>2 + 

O 

[1 

CO 

u* 

(21c) 

m 2 ' - V, - P| 

S' + “2 

) * t 2 = 0 

(2 Id) 

M 3 , + V 2 + P ( 

1 > 
y 2 -"3, 

) + t 3 = 0 

(21e) 


with boundary conditions at the ends of the beam being 


XJ 1 = 0 

or 

N - P - N = 0 


to 

I! 

’ O 

or 

nu 2 ' + v 2 - pc 3 

- v 2 -° 

CO 

II 

o 

or 

NU 3 ’ + V 3 + Po> 2 

O 

11 

CO 

!> 

i 

o 

li 

CM 

3 

or 

O 

11 

CM 

i§ 

1 

CM 

2 


o 

1! 

CO 

3 

or 

O 

11 

CO 

12 

CO 

2 



( 22 ) 


The assumption of zero twist allows the underscored terms in equations (13) to be 
dropped. Since the first of equations (13) is still nonlinear, equations (21) remain, in one 
sense, nonlinear, because they retain their ability to handle problems involving large 
lateral displacements. However, whether or not the first of equations (13) is linearized, 
equations (21) are "formally" linear because equation (21a) may be solved without regard 
for the lateral loads and deformations. Equations (21), which will be referred to as 
linear, are equivalent to the well-known Timoshenko beam equations (see ref. 8), with the 
additional effects of axial force, internal pressure, and large lateral displacements 
accounted for. 


ILLUSTRATIVE EXAMPLES 

Two sample problems have been solved to illustrate the usefulness of the simplified 
system of equations. The first problem is that of the buckling of a pin-end column under 
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a compressive axial force. The second problem is concerned with finding the displace- 
ments and rotations of a pin-end beam column under a compressive axial force and a 
linearly varying transverse load. For simplicity, displacement and rotation only in the 
XY-plane are considered. 

Problem 1.- The configuration to be analyzed is shown in figure 2. 



If, in equations (21), 


qi . q 2 = q 3 = T 2 = T 3 = m 2 = m 3 = co 2 = U 3 = ° 


and 


then, from equation (21a), 


N = -F 


N = P - F 


and equations (21b) and (21e) become 

(P + Gurt - F)U 2 m - (P + G7rrt)o> 3 ' = 0 

Emr 3 ta> 3 " - (P + G7rrt)(u 2 ' - = 0 

The boundary conditions are 

U 2 (0) = U 2 (L) = w 3 ’(0) = o> 3 '(L) = 0 
Integration of equation (23a) and substitution into equation (23b) give 

ct > 3 +Po> 3 - £_ 


(23a) 

(23b) 


(24) 


(25) 
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where 


^2 _ F / P + G7rrt \ 
ETrr^typ + Gwrt - F / 


and C is a constant of integration. Then 

C 

cog = Bj sin /3x + Bg cos )3x + — 

Use of the boundary condition <Ug'(0) = ci»g'(L) = 0 (eq. (24)) yields 


( 26 ) 


and, for a nontrivial solution 


B. = 0 


/3L = 7T 


(27) 


from which equation the lowest critical value of F is found. Substitution for /3 into 
equation (27) and substitution of I for 7 rr^t gives 


F = 
cr 


El ^P + Gvrrt) 


El Z_ + p + Girrt 


(28) 


This result is in agreement with the results presented in reference 9. 

In equation (28) may be seen the relationship which obtains between the bending 
stiffness and the shear stiffness (including shear stiffness provided by internal pressure) 

in the expression for the critical axial compressive force. Denoting the Euler buckling 

„2 

load El — by Kg and the shear stiffness P + G 7 rrt by P + Kg transforms equa- 
L 

tion (28) to 


F = 
cr 


k b( p + K s) _ 

K B + P + K s 


K 


B 


P + K c 


1 + 


P + K c 


1 + 


P + K c 


K 


B 


(29) 
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the shear buckling load. Equation (29) shows that internal pressure increases the 
resistance to transverse shear deformation. In fact, for inflated bodies made of single- 
layer fabrics for which a low effective shear modulus may be expected, the pressure con- 
tribution may well predominate. (See ref. 10.) 

Problem 2 .- The second example is an application of the present results to a 
simplified paraglider problem. It was shown in reference 5 that the stresses normal to 
the leading edges of the paraglider sail vary linearly along the edges. Then if the 
leading-edge boom is considered to be a long inflatable cylindrical beam, and the pay- 
load Q is positioned directly under the centroid of the linearly varying load on the 
boom by two shroud lines of lengths l ^ and l ^ (see fig. 3), the problem may be 

defined in the following manner: The condi- 

tions of problem 1 still hold, except that 



where, from considerations of equilibrium 
and the geometry of the boom-shroud- 
payload configuration, 




Figure 3.- Pin-end beam column under compressive axial force 
and transverse load. 
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and 


h-\h 2 4) 


1/2 


with l j given ^ > | l| . 

Equations (21b) and (21e) become 


(P + K S -N)U 2 ”-(P + K S ) W 3' = 


_ 2Qx 
T 2 


and 


EI" 3 " + ( P + K S )(U 2 ' - « 3 ) = 0 

with the boundary conditions 

U 2 (0) = U 2 (L) = w 3 ’(0) = u> 3 '(L) = 0 

Solution of equations (30) and (31) subject to the boundary conditions yields 


( 30 ) 


(31) 


(32) 


and 


U 


2EIQ/sin Ax 
N 2 L \sin XL 



+ 


Qx/x 2 

3N\L 2 


(33) 




= ^ 
N 


COS Xx | l/x 2 

- 2 - A 

XL sin XL 2 \L 2 

x 2 l 2 sj 


(34) 


where 

.2 S (P + K^ 

EI^P + Kg-Nj 

Expressions for the shear V 2 and the moment M 3 may be obtained by use of equa- 
tions (33) and (34) along with equations (13) and (15): 


V 2 = 


2QEI 


NL 2 


"ir 


X L cos Xx 
sin XL 


- 1 


(35) 
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Another quantity of interest is the relative axial displacement of the ends of the 
cylindrical beam. From equations (13a) and (15a), if N = P - N, 

U ’ - P -N 1 tj - 2 

1 EA 2 U 2 

where A = 27;rt. With the condition that U^(0) = 0, the axial displacement is 


U i = 


p - n x i r x rr > 2 


EA 


i r 

"2J 0 


U 2 dx 


(37) 


Substitution for U 2 (x) and integration yield for the relative axial displacement of the 
ends 


U X (L) = P - N - L - L< 


EA 


+ 2 1 - 


2S?. + 2 


|45N 2 


/ EIQ 
\N 2 L 2 


A 2 + it cot XL - 1 


2 sin^XL 2 


N 


f x 2 l 2 ^ 

XL cot XL + - 1 


P + K c 


(38) 


If the transverse shear stiffness is taken to be infinite, these results can be shown to be 
in agreement with the results for the classical beam column. 


CONCLUDING REMARKS 


Under the assumption that the deformations of long pressurized cylindrical mem- 
brane beams can be accurately expressed in terms of translations and rotations of a 
cross section, nonlinear equilibrium equations have been derived for the bending, twisting, 
and stretching of such beams. It is found that for problems in which twist about the lon- 
gitudinal axis can be neglected the equilibrium equations are linear. (Twisting will not 
occur in beams of circular cross section when they are loaded as a beam and/or column.) 
These linear equations are seen to be equivalent to the Timoshenko beam equations with 
the transverse shear stiffness augmented by the appropriate internal pressure contribu- 
tion and with axial force accounted for. 

Two illustrative problems are analyzed by using the linear equations. The first 
problem is that of the buckling of a simply supported pressurized column under axial 



compression. The results are seen to be in agreement with results obtained by another 
method. 

The second problem, typical of those which could arise in the use of pressurized 
membrane beams as booms for a paraglider or similar vehicle, is concerned with the 
combined bending and axial compression of a simply supported inflated membrane 
column. Expressions are obtained for lateral deflection, cross-section rotation, shear 
force, bending moment, and column shortening. If the transverse shear stiffness of the 
inflated beam column is assumed to be infinite, the results can be shown to agree with 
the results for the classical beam column. 

The analysis presented herein is valid for slender beams having significant wall 
bending stiffness, with or without internal pressure, provided local buckling of the wall 
does not occur. Since the theory contains no provision for the effects of local buckling, 
it is applicable to membrane cylinders only when the internal pressure and applied 
forces are such that the cylinder wall is in tension. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., February 14, 1966. 
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